1 



m 
o 
o 

(N 
-i— > 
O 

O 

oo 



> ■ 
<n : 

in ■ 
o ■ 

m : 
o ■ 

^ ' 

6 
& 

in 
c3 
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Jacob Bekenstein's identification of black hole event hori- 
zon area with entropy proved to be a landmark in theoretical 
physics. In this paper we trace the subsequent development 
of the resulting generalized second law of thermodynamics 
(GSL), especially its extension to incorporate cosmological 
event horizons. In spite of the fact that cosmological hori- 
zons do not generally have well-defined thermal properties, 
we find that the GSL is satisfied for a wide range of models. 
We explore in particular the case of an asymptotically de Sit- 
ter universe filled with a gas of small black holes as a means 
of casting light on the relative entropic 'worth' of black hole 
versus cosmological horizon area. We present some numeri- 
cal solutions of the generalized total entropy as a function of 
time for certain cosmological models, in all cases confirming 
the validity of the GSL. 



I. INTRODUCTION 

A key development in the history of physics came with 
Jacob Bekenstein's identification (Bekenstein, of en- 
tropy Shh with black hole event horizon area Abh, 



Sbh oc k B A 



bh 



(1) 



where k B is Boltzmann's constant. That there was a 
close analogy between black hole event horizons and en- 
tropy was implied by Hawking's area theorem, which 
states that subject to certain reasonable physical condi- 
tions (most notably that energy cannot be negative) the 
total horizon area is a non-decreasing function of time 
(Hawking, 0). However, if a black hole were totally 
black it would have a zero temperature. Assuming the 
general relationship: entropy = energy/temperature, it 
would seem that the entropy of a black hole diverges. A 
similar conclusion follows from information theory. If the 
black hole forms from the implosion of a ball of matter, 
the information lost behind the event horizon is roughly 



N bits, where N is the number of particles in the ball. 
As classical physics imposes no lower bound on the mass 
of a particle, there is no upper bound on N. Then iden- 
tifying entropy with information loss confirms that Sbh 
diverges. 

It was Bekenstein's suggestion, with the encourage- 
ment of John Wheeler, that quantum mechanics re- 
moves the divergence by placing a lower bound on the 
mass of the particles that go to make up the black hole 
(Bekenstein, pjj). In order to confine a particle to the 
Schwarzschild radius, its Compton wavelength must be 
less than the size of the hole, from which Bekenstein con- 
cluded that the constant of proportionality in Eq. ^ in- 
cludes the factor l/fr. These essential ideas were later 
confirmed following the application of quantum field the- 
ory to the formation of black holes by Hawking , who 
discovered that a Schwarzschild black hole of mass Mbh 
and surface gravity Kbh radiates with the temperature 



Tbh = h/(2Trk B c)K hh 

he 3 

~ 8Trk B GM hh 
This fixes the constant of proportionality in Eq. ^ 



(2) 
(3) 



Sbh 



A 



k B c 3 
\A hh 



bh 



(4) 
(5) 



where Eq. [3] uses units with h — c=G = k B = l. We 
shall adopt these units henceforth. 

Following these developments, the second law of ther- 
modynamics could then be generalized to include cases 
where black holes exchange heat and energy with their 
environment (Hawking, 



S 



bh 



s m > 



(6) 



where S m is the entropy of the matter and an overdot rep- 
resents differentiation with respect to proper time. The 
fact that black holes radiate implies that they can lose en- 
ergy and shrink, in violation of Hawking's area theorem. 
The strictures of the theorem are evaded because the 
quantum state permits a flux of negative energy into the 
hole (Davies, Fulling, and Unruh, However, the ther- 
mal radiation emitted by the black hole always raises the 
entropy of the environment by at least as much as the loss 
of horizon entropy caused by the shrinkage. The general- 
ized second law (GSL) was extended to black holes with 
rotation and electric charge in a straightforward manner 
(Hawking, 0). 
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FIG. 1: Scale factor is plotted against time for a flat, ra- 
diation filled, FRW universe fEq. If 311 . Hubble's constant is 
taken to be Ho = 70kms~ 1 Mpc~ 1 , and the cosmological pa- 
rameters Qm = 8irGpo/3Ho and Qa = A/3Hq are given their 
observationally favoured values 0.3 and 0.7 respectively. The 
expansion initially decelerates while the gravity of the radia- 
tion dominates, then begins to accelerate as the cosmological 
constant takes over. In this universe, as in all eternally ex- 
panding FRW universes with non-zero A, the expansion tends 
toward exponential at late times, a(t) oc e Ht . 



II. EXTENDING THE GENERALISED SECOND 
LAW OF THERMODYNAMICS TO 
COSMOLOGICAL HORIZONS 



Shortly after the establishment of black hole horizon 
entropy, it was proposed that the GSL be extended to in- 
clude the event horizon area of de Sitter space (Gibbons 
& Hawking, 6]). This cosmological model has attracted 
interest in the last twenty years because, if the inflation- 
ary scenario is correct, the very early universe experi- 
enced a period of exponential expansion that approxi- 
mated a portion of de Sitter space. Moreover, recent cos- 
mological observations suggest that the cosmological con- 
stant A is non-zero and positive (see Lineweaver and 
Sievers et al. @ for recent reviews) — all ever-expanding 
Friedmann-Robertson- Walker (FRW) universes of this 
class tend towards de Sitter space at late times. Hence 
universes with event horizons have recently experienced 
a resurgence of interest. 

De Sitter space has a time-independent event horizon 
at the proper distance r c = 1/H = \/3/A (subscript c 
stands for cosmological event horizon). Given Hubble's 
law v — Hr, a simple rearrangement shows that the event 
horizon is at the distance at which, using the above def- 
inition of v, comoving test particles recede at the speed 
of light. In a de Sitter universe particles receding faster 
than light are beyond our view (but not in general for 
FRW universes (Kiang, ||)). 

The thermal properties of de Sitter space may be de- 
duced from the fact that the vacuum Green function 
for conformally-invariant fields propagating in a de Sit- 



ter background spacetime are periodic in imaginary time 
with a periodicity, 2tt/k c , corresponding to a thermal 
state with temperature 



T, 



dcS 



K C /27T, 

-(- 

2ir I 3 



1/2 



(7) 
(8) 



where re c = (A/3) 1 / 2 is the surface gravity of the cos- 
mological event horizon. The same result may be estab- 
lished by considering the response of a particle detector 
travelling along a geodesic in a de Sitter-invariant vac- 
uum state; the detector behaves as if immersed in a bath 
of thermal radiation with temperature given by Eq. 
Gibbons and Hawking [|| therefore made the reasonable 
suggestion that the GSL be extended to de Sitter space. 
(For more general derivations of horizon entropy see the 
recent pioneering work of Padmanabhan [23l 123. Hj|.) 
There were, however, some considerations that made the 
identification of de Sitter event horizon area with entropy 
less compelling than in the black hole case. First, de Sit- 
ter horizons do not radiate in the manner of black holes. 
In spite of the thermal response of a particle detector, the 
expectation value of the stress-energy momentum tensor 
of the de Sitter-invariant massless scalar field vacuum 
state is not that of a bath of thermal radiation. It is 
given instead by (see, for example, Birrell and Davies, 

M) 

A4 

T ^ = 8640^^- (9) 
This corresponds to the stress-energy-momentum tensor 
associated with a cosmological constant, and so merely 
renormalizes A. Second, in the black hole case the ob- 
server lies outside the horizon, whereas in the cosmo- 
logical case the horizon not only envelops the observer, 
its location is observer-dependent. There is therefore 
no apparent source of the radiation; it appears homo- 
geneous. Third, there is no parameter in de Sitter space 
that corresponds to the mass of the black hole, making 
it problematic to balance the books in a trade-off of en- 
ergy and entropy between the de Sitter horizon and any 
matter present. Finally, a key heuristic in Bekenstein's 
derivation of Eq. ^ was the association of information 
with entropy. The black hole irreversibly swallows the 
information of the body that implodes to form it, and 
the entropy of the horizon may be thought of as a mea- 
sure of this lost information. In the case of de Sitter 
space, the region beyond the horizon is infinite, and can 
accommodate an infinite quantity of information. 

In spite of these shortcomings, we shall show that the 
GSL can be consistently applied to de Sitter horizons 
that exchange heat and energy with both other horizons 
and with matter. 

To investigate the entropy and energy budgets when 
the de Sitter horizon exchanges heat with other systems, 
we first consider the case of a FRW universe with a posi- 
tive A containing uniform thermal radiation with temper- 
ature T > Td s = H/2it (which corresponds to a typical 
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wavelength of the radiation being less than the horizon 
radius). Consider a FRW universe model with scale fac- 
tor a, and Robertson- Walker metric 

ds 2 = dt 2 - a 2 (t) [d X 2 + S k {x) 2 (d6 2 + sin 2 6d<p 2 )] (10) 

where Sk(x) — sinhx, x, si n X f° r & = —1,0,1 respec- 
tively. The Friedmann equations may be written 



P 

3H 2 



-3.ffOH 
8np + A 



P) 

-3k/a 2 



(11) 
(12) 



where p and p are the density and pressure of the cos- 
mological fluid respectively, and we take p = p/3 for 
radiation. In the spatially flat case, k = 0, Eqs. El and 
E| may readily be solved to give the scale factor 



a(t) = 



\ A 



1/4 



sinh 




(13) 



where po is the present day radiation density and p — 
poa~ 4 . Eq. El is plotted in Fig.Q] from which it can be 
seen that the universe starts out with a standard hot big 
bang, and approaches de Sitter space as t — > oo. 

Consider a small departure from de Sitter space, such 
that the distance to the cosmological horizon, r c , is ap- 
proximately 1/H, corresponding to an entropy of S c ~ 
ir/H 2 . This will be a good assumption at late times. The 
rate of change of cosmological horizon entropy is then ap- 
proximately 



S c w -2irH/H 3 = 2tt(16ttp)/3H 3 



(14) 



using Eq. 1121 This rise in horizon entropy will be 
traded against a loss of radiation entropy as the radiation 
streams away across the de Sitter horizon. The radiation 
energy density can be expressed as p = crT 4 where the 
radiation constant a = 7r 2 /15 in these units. The entropy 
density is s — (4/3)/oT _1 , so the total radiation entropy 
within a horizon volume is given by 



(15) 



Using Eq. (11), the rate of loss of radiation entropy is 



^^(S/ip-^pH- 3 - 3pV 4 HH- 



-IQttu 1 ' 4 



p 3 / 4 H- 2 



l ) (16) 
(17) 



if H is small. 

satisfied is 



The condition for the second law to be 



S c + S x > 0. 



(18) 



Comparing Eq. El with Eq. El one sees that the gain in 
horizon entropy exceeds the loss of radiation entropy if 
p x l A > H/2n. This is just the condition that the radi- 
ation temperature be greater than the horizon tempera- 
ture, as assumed. 



In the foregoing we have assumed a small departure 
from de Sitter space. Unfortunately when we relax this 
condition the distance to the event horizon cannot be 
solved exactly, but we show a numerical solution in Fig. [21 
From this it may be seen that the horizon entropy always 
rises faster than the loss of radiation entropy, confirming 
the validity of the GSL even at early times. In arriving 
at these results we define the (time-dependent) horizon 
radius as 



a(t) 



dt' 

WY 



(19) 



It is possible to prove an analogue of the Hawking area 
theorem for cosmological horizons (Davies, ). Re- 
call that it was this theorem that originally motivated 
the association of entropy with black hole horizon area. 
Consider the class of homogeneous isotropic models filled 
with a fluid of pressure p and energy density p and met- 
ric given by Eq. 1101 It is convenient to work with the 
conformal time parameter defined by 



V = 



dt' 



(20) 



If \rj\ < go for a particular model universe, then that 
universe possesses an event horizon. 

Comparing Eq. 1191 and Eq.JSOIwe see that the distance 
to the cosmological event horizon is given by r c = —at]. 
When the area of the event horizon is calculated in the 
k = ±1 cases we must take into account curvature terms. 
We prefer to work with the corrected distance 



-aS k (rj)- 



(21) 



With this definition, the area of the event horizon always 
increases if 



> 



(22) 



where a prime indicates differentiation with respect to 
r\. Performing this differentiation on Eq. |^ gives the 
condition 



(23) 



where Tk( X ) = tanh.%, x,tan% for k = —1,0,1 respec- 
tively. It is interesting to give a physical interpretation 
of Eq. El For the flat (k = 0) case it implies that the 
event horizon area does not decrease as long as d > —Iff], 
or r c > 1/H. That is, the event horizon area always in- 
creases as long as it is farther away than the Hubble 
sphere (the distance at which the recession velocity, us- 
ing the definition v — Hr, is equal to the speed of light). 
Note that the inequality allows the event horizon to be 
more distant than the Hubble sphere: in most FRW uni- 
verses we can observe objects which are receding faster 
than the speed of light. The condition is slightly varied 
for the curved cases, because in these cases it is some- 
times possible for objects receding less than the speed of 
light to be beyond the event horizon. 
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We now assume that the cosmological fluid is subject 
to the dominant energy condition p + p > 0. This is the 
same condition as used in the black hole area theorem 
— the one that is violated by quantum effects, allowing 
black holes to evaporate. It then follows from Eq. ^] 
that p < 0, from which inequality Eq. |23 readily follows 
(Davies, [n)). 

In the cosmological case, there is no direct analogue of 
the evaporation of a black hole. However, one may relax 
the dominant energy condition by considering cosmolog- 
ical fluids with pressure p = (7 — l)p and bulk viscosity 
up > 0. It is then possible (Davies, [l2]) to have a com- 
bined effective pressure of the fluid plus the cosmological 
contant in excess of minus the combined energy density. 
Equation 1141 is then modified to 



fall at the roots of the quartic equation 



S c = 8ir 2 (j -3Ha)p/H 3 



(24) 



from which it follows that the horizon entropy will de- 
crease if 7 < 3Ha. This does not, however, necessarily 
signal a violation of the GSL. The presence of viscosity 
implies the generation of heat entropy. A simple calcu- 
lation (Davies, shows that the fluid generates heat 
entropy at a rate 



S n 



9H 2 apa 3 
T 



(25) 



where T is the temperature of the fluid and we identify 
a 3 at an instant of cosmic time, t, with a horizon vol- 
ume. Comparing Eq. [31 and Eq. [25] for the case 7 = 
(corresponding to minimal departure from the condition 
p + p = 0) shows that the total entropy is unchanged for 
the equilibrium case 



T 



H 

2^ 



T, 



dcS 



(26) 



thus respecting the GSL. The case T > T^eS requires a 
more complicated treatment (Davies, [Lql but also con- 
forms with the GSL. Several exact solutions of cosmolog- 
ical models with viscosity have been discussed in detail 
by Barrow and Davies [T^|. 



III. BLACK HOLE-DE SITTER SPACETIMES 

Special interest attaches to cosmological models that 
contain both black holes and cosmological horizons. As 
these spacetimes evolve, there may be an exchange of 
horizon area between the two. The question then arises of 
whether cosmological horizon area has the same 'worth' 
as black hole horizon area. If the loss of black hole hori- 
zon area was not compensated by at least as great gain 
in cosmological horizon area, for example, then the GSL 
would be threatened. 

A general class of black- hole-de Sitter solutions includ- 
ing mass M, charge Q, and angular momentum J was 
provided by Carter |15| . The horizons for these solutions 



- J 2 I r 2 



2Mr 



j 



A 



= 0. (27) 



If the four roots are given by r ++ > r + > r_ > r 

then r ++ is identified with the cosmological event hori- 
zon, and r+ with the outer horizon of the black hole. For 
an observer in the region outside the hole, but within the 
cosmological event horizon, the total horizon area is 



A tot = 4tt(t 



2J 2 ) 



(28) 



(Gibbons and Hawking, 6]). This spacetime is thermo- 
dynamically unstable: the black hole temperature is al- 
ways greater than the cosmological horizon temperature. 
The black hole slowly evaporates, and the resulting ra- 
diation eventually passes over the cosmological horizon. 
At late times the solution settles down to de Sitter, with 
horizon area 



Aics - 127TA- 1 . 



(29) 



It follows from the algebraic conditions on the four roots 
of the quartic (Eq. E3 that 



(30) 



from which it readily follows that Eq. [35] is greater than 
Eq. EHl in conformity with the GSL (Davies, flfil|). For 
black holes with charge Q and angular momentum J the 
proof is limited to the parameter range A" 1 > 3M 2 > 
3J 2 . This is roughly the range for which the cosmological 
horizon lies outside the black-hole horizon and the central 
mass is not a naked singularity. 

Davies, Ford and Page [17| extended the foregoing 
'before-and-after' analysis by considering the continuous 
exchange of small quantities of heat between the black 
hole and the cosmological horizon. They introduced 
a further generalization by considering a massive shell 
around the black hole. This has the effect of depress- 
ing the temperature of the hole. A sufficiently massive 
shell lying sufficiently close to the black hole horizon can 
produce a backflow of heat from the cosmological hori- 
zon into the black hole. It turns out that even in these 
circumstances the GSL is satisfied. 

As a final refinement of the investigation of the trade- 
off between black hole and cosmological horizons, we con- 
sider a FRW model universe with positive cosmological 
constant and a fluid consisting of uniform dust composed 
entirely of identical slow-moving small Schwarzschild 
black holes of mass Mbh and number density rtbh- We 
ignore the Hawking effect (it could be suppressed if nec- 
essary by considering maximally charged black holes) and 
the small correction to the black hole radius implied by a 
non-zero A. In this case the Friedmann equations, Eg. 1111 
and Eq. ^] can be written 



Pbh 

m 2 



-3Hp hh 
87rpbh + A 



3k/a 2 



(31) 
(32) 



500 



,-400 



300 



. 200 



100 



III 

(n M ,n A ) = (0.3,0.7) 




- 

— Cosmological entropy 




— Radiation entropy 




— Total entropy 



10 

Cosmic time, 



(Gyr) 



15 



20 



500 



400 



co 300 



100 



("„." A ) = (0.3,1.4) 




r bh =1.0 cGyr 

/ / 




- 




— Cosmological entropy 




Black hole entropy 




— Total entropy 



10 15 
Cosmic time, t, 



20 



(Gyr) 



25 



30 



FIG. 2: The entropy of the cosmological horizon, the entropy 
of the radiation within the horizon and the combined total 
entropy are plotted against time for a radiation dominated 
FRW model with (Qm.^a) = (0.3,0.7). At all times the 
loss of radiation entropy across the cosmological horizon is 
less than the increase in cosmological horizon entropy, so the 
GSL is satisfied. 



FIG. 3: Entropy variation as black holes cross the cosmologi- 
cal event horizon. The model shown is a FRW universe with 
(Hm, Oa) = (0.3, 1.4) in which the matter density consists en- 
tirely of a comoving gas of black holes. This is an example of 
a closed universe, and the cosmological event horizon appears 
out of the antipode at a finite time — prior to this the black 
hole entropy is constant and the cosmological horizon entropy 
is zero. 



where p hh = M hh n bh . 

Restricting for the moment to small perturbations 
about de Sitter space, we may proceed in the same fash- 
ion as the calculation in Eqs. 1141181 above. Eq. ^] gives 
the rate of change of the cosmological horizon entropy, 
S c = A c /A. The total black hole area within the cosmo- 
logical horizon is, 



A b ° h tal = A hh n hh V c 



64ir 2 M hhPhh 
3iJ 3 

The rate of change of black hole area is therefore, 



.4 



\ total 
bh 



647r 2 M bh/ o bh 
H 2 



(33) 
(34) 

(35) 



where we have used Eas. 1311 and 1321 and once again ig- 
nored the H term. For the GSL to hold we need, 



i b ° h tal + Ac > 0. 



Using Eas. 1141 and 1351 this inequality becomes 



2Af bh 



1/H 



(36) 



(37) 
(38) 



Thus the GSL holds as long as the black holes are smaller 
than the cosmological event horizon. 

We have been unable to solve inequality 1201 exactly for 
all FRW universes, but we give a sample numerical solu- 
tion for large departures from de Sitter space in Fig. 
Note that the total entropy is a monotonically increasing 
function of cosmic time t at all times. We shall report 
elsewhere (Davis et al. [3) on a broader range of nu- 
merical solutions of black hole cosmological models. 



IV. RELATION TO OTHER WORK 



In this review we have concentrated on several explicit 
examples, including some numerical solutions, of exten- 
sions of Bekenstein's generalized second law of thermo- 
dynamics to cosmological horizons. This work comple- 
ments some recent theorems that prove more general but 
less explicit results that have a bearing on the GSL. For 
example, for black hole-de Sitter spacetimes, Shiromizu 
et al. .19] show that the black hole event horizon area is 
non-decreasing in asymptotically de Sitter space times, 
while Hayward et al. |2J| show that the black hole event 
horizon area is bounded by 47r/A. 

The cosmological horizon area in models with a 
positive cosmological constant has been considered by 
Boucher et al. |2lJ who show that the horizon area is 
bounded by 127r/A on a regular time-symmetric hyper- 
surface, while Shiromizu et al. [Tg| show that the hori- 
zon area is bounded by 127r/A on a maximal hypersur- 
face. Neither of these proves a bound in a non-stationary 
asymptotically de Sitter universe such as a FRW uni- 
verse. 

Maeda et al [2^ extend the work of Davies ^l| by 
showing that the cosmological event horizon area doesn't 
decrease in any asymptotically de Sitter spacetimes. 
They also show that the de Sitter horizon is the up- 
per limit of horizon size for any cosmological model with 
nonzero A. Our results both illustrate these theorems and 
demonstrate that, for certain specific models, the GSL is 
satisfied not just asymptotially, but at all times. 
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Postscript. 

The reference to Padmanabhan [23l l24l |25| was inad- 
vertantly omitted from the published version. 
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